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The transition temperature for a dilute, homogeneous, three-dimensional Bose 
^ ■ gas has the expansion T^ = ro{l+cian^/^ + [c2 ln(an^/^)+C2]a^n^/'^+0(a^n)}, 

'^ . where a is the scattering length, n the number density, and Tg the ideal gas 

d , result. The first-order coefficient ci depends on non-perturbative physics. In 

this paper, we show that the coefficient c'2 can be computed perturbatively. 
,-^ ' We also show that the remaining second-order coefficient C2 depends on non- 

|Z| ■ perturbative physics but can be related, by a perturbative calculation, to 

quantities that have previously been measured using lattice simulations of 

three-dimensional 0(2) scalar field theory. Making use of those simulation 
(NJ ; results, we find T^ ~ Tojl + (1.32 ± 0.02) an"'^ + [19.7518 hi{an^/^) + (75.7 ± 

>■ 0.4)]a2n2/3 + O(a3n)}. 
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!>. ■ I. INTRODUCTION 

o 

Long-distance physics at a second-order phase transition is generically non-perturbative. 
For this reason, researchers have found it non-trivial to compute corrections to the ideal 
gas result for the critical temperature T^ for Bose-Einstein condensation (BEG) of a dilute, 
homogeneous Bose gas in three dimensions. It is currently understood that the correction 
^ . ATc = Tc — To to the ideal gas result Tq behaves parametrically as 

O ■ ATc 1/3 



^0 



cian''^ (1.1) 



in the dilute (or, equivalently, weak-interaction) limit, where ci is a numerical constant, 
5^ I a is the scattering length, which parameterizes the low energy 2-particIe scattering cross- 

section, and n is the density of the homogeneous gas. We will assume that the interactions 



are repulsive (a > 0). A clean argument for ( |1 . 1|) may be found in Ref. ||l|, which also 
shows how the problem of calculating the constant Ci can be reduced to a problem in three- 
dimensional 0(2) field theory. Recent numerical simulations of that theory have obtained 
the results Ci = 1.29 ± 0.05 § and Ci = 1.32 ± 0.02 ||,|. 

In this paper, we shall extend the result for Tc{n) to second order in a for a homogeneous 
Bose gas. This is also the relationship between T^. and the central number density for a 



Bose gas in an arbitrarily wide trap. (In contrast, the relationship between Tc and the total 
number A^p of particles in a trap depends on somewhat different physics. A second-order 
result for Tc{Np) in an arbitrarily wide trap may be found in Ref. [^].) 

In the homogeneous case, Holzmann, Baym, and Laloe ^ have recently argued that a 
logarithmic term appears at second order, 

AT 

f^ ^ CianV3 + c',aW^ Inian'^') + ■ ■ ■ , (1.2) 

and they made a rough estimate of the coefficient Cg using large- A^ arguments. A similar 
logarithm has been found for Tc{Np) in the case of trapped gases 0. We will show that, in 
contrast to Ci, the coefficient c'a of the logarithm can be computed exactly using perturbation 
theory. Our result is 

_ 64H(^) ,j 3, 

- 'c(i)r" ' ' ' 

where ({s) is the Riemann zeta function. We will compute this result by performing a second- 
order perturbative calculation to match the physics of the transition onto that of three- 
dimensional 0(2) scalar field theory. The same matching calculation will also perturbatively 
determine the relationship between the non-logarithmic term at second order and certain 
non-perturbative quantities in 0(2) scalar theory which have been previously measured in 
lattice simulations. As a result, we will determine all the coefficients in the second-order 



expansion 

AT, 



To 



cxan^l'^ -F 



,2^2/3 



^ ln(ani/3) + 4' aV/=^ + ©(a^n) , (1.4) 



[where the notation 0[a^n) is not intended to make any particular claim about what powers 
of logarithms might appear at third order] . 

We should emphasize that when we refer to the "first order" and "second order" terms 
in ( |1.4| ), we do not mean first and second order in perturbation theory. Perturbation theory 
breaks down for these quantities, and that breakdown manifests as the appearance of infrared 
(IR) infinities beyond a certain order in the perturbative expansion. 

A portion of the required perturbative matching calculations has already been performed 
in Ref. |]^, which also gives a discussion of the philosophy and methods of perturbative 
matching calculations between Bose gases and three-dimensional 0(2) field theory at the 
phase transition. However, we find it convenient to use slightly different conventions than 
Ref. . For the sake of introducing conventions and notation, and for the sake of making 
this article somewhat self-contained, we will use the remainder of this introduction to briefly 
review matching and the different distance scales associated with physics at the transition. 
In section |I|, we will flx the ultraviolet (UV) regularization and renormalization schemes we 
will use for our calculation. We then proceed to do the matching calculations and assemble 



all the matching results in section |ITI| , though we leave the details of the more intricate 
diagrammatic calculations for later sections. In section ^, we then put together our flnal 
results for the relationship between T and n at the transition. Section gives the details of 



how to calculate the most complicated diagram that was needed for matching. Section VI 



reproduces some previous results related to the critical value fJ,c{T) of the chemical potential 
01 in a new form that is needed for our analysis. Finally, section [VII| explains how our result 
for the second-order logarithm is modified in theories with A^ fields, for the sake of readers 
who may wish to compare exact results to the approximate large- A^ analysis of Holzmann, 
Baym, and Laloe ^ . A very brief outline of how to calculate a few simple finite-temperature 
integrals in dimensional regularization is left for an appendix. 

For the calculations in this paper, it will be convenient to follow Baym et al. |]I| and cal- 
culate the critical density nc{T) as a function of T rather than the critical temperature Tc(ra) 
as a function of n. We will then obtain the formula for Tc{n) by inverting the relationship. 
The ideal gas result no{T) for ric is 

C(-) 
no{T) = ^, (1.5) 

where 



is the thermal wavelength. (In this paper, we work in units where h = k^ = 1, where /cb 
is Boltzmann's constant.) The diluteness condition an^^^ <^ 1 for the expansions discussed 
above can therefore alternatively be expressed as 

a<A (1.7) 

at the transition. 

A. Overview of matching to 3-dimensional 0(2) theory 

Baym et al. [|l| were the first to use an effective three-dimensional 0(2) scalar field theory 
to study non-universal long-distance physics of the BEC transition of a dilute Bose gas. A 
more systematic discussion of how to match the parameters of the 0(2) theory to the original 
problem, in order to study interaction effects beyond leading order, may be found in Ref . . 
Here, we will briefly review these issues in preparation for doing the matching calculations 
that we will need to obtain Tc{n) at second order. Perturbative matching calculations, which 
allow effective theories to be used to calculate non-universal quantities, can be performed 
whenever the s/iort-distance physics described by the effective theory is perturbative, at a 
scale where the effective theory is still applicable. Such calculations have a long history 
that includes lattice field theory [0, Bose condensation at zero temperature 0, relativistic 
corrections to non-relativistic QED [^J, heavy quark physics [|r^, ultra- relativistic plasmas 
TT| , and non-relativistic plasma physics [jl2[. For a general discussion, see also Ref. [^ . 



The starting point is the well-known description of a dilute Bose gas by a second- 
quantized Schrodinger equation, together with a chemical potential /i that couples to particle 
number density 4'*'^^ ^ind a IV'I*^ contact interaction that reproduces low-energy scattering 



I^ . The corresponding Lagrangian is 

1 _o \ , 27ra 



C = r(^^t + 7^V' + ^^)^|J--^ (r^)'. (1.8) 

V 2m J m 



This effective description is valid for distance scales large compared to the scattering length 
a. Corrections to this description, due to the energy dependence of the cross-section or 3- 
body interactions or so forth, do not affect Tc{n) at second order (see section |IV|) . At finite 
temperature, it is convenient to study ( |1.8|) using the imaginary time formalism, in which t 
becomes —it and imaginary time r is periodic with period (3 = 1/T. The imaginary-time 
action is then 



Si= dr I d^x 



2m J m 



:i.9) 



We shall call this the 3+1 dimensional theory, referring to three spatial dimensions plus one 
(imaginary) time dimension. The expectation value of the number density is given by 



n = {1^*1^). 



;i.io) 



As usual, the field ip can be decomposed into imaginary-time frequency modes i/jj with 
discrete Matsubara frequencies ujj = 27ij/P, where j is an integer. If we ignore interactions 
for a moment, and treat /i as small, then non-zero Matsubara frequency modes in ( p..9|) 
are associated with a correlation length of order {2miUj)~^^'^ ^ A, where A is the thermal 
wavelength (|1.6| ). Near the transition, at distance scales large compared to the thermal 
wavelength A, all the modes with non-zero Matsubara frequencies decouple, leaving behind 
an effective theory of just the zero- frequency modes V'o(x). If one were to naively throw 
away the non-zero frequency modes from the original 3+1 dimensional action ( |1.9D , it would 
reduce to 



Si^/3 d^x 






/U hAo 



2-Ka 
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(^o>o) 



1.11) 



This is the rough form of the effective 3-dimensional theory of ipo, which describes long- 
distance physics at the transition for time-independent quantities. However, completely 
ignoring the effects of non-zero frequency modes was an oversimplification. In field theories, 
the short- distance and/or high-frequency modes do have effects on long-distance physics, but 
those effects can be absorbed into (1) a modification of the strengths of relevant interactions 
(in the sense of the renormalization group) between the long-distance/zero-frequency fields, 
and (2) the appearance of additional marginal and irrelevant interactions between the long- 
distance/zero-frequency modes. The latter effect will not be relevant at second order for 
Tc{n) (see section fV|). Because of the first effect, the correct three-dimensional effective 
theory is of the more general form 



S3 = P d^x 



ro[- 



2m 



V" - /i3 hAo + Z, 



27ra 



m, 



iro^of + h 
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where the difference of the (^o-independent) parameters Z^, /is, Za, and /s from the naive 



values Z 



V> 



/^3 



/i, and /s = of (|1.11|) incorporates the effects of short- 



distance physics on long-distance physics. These parameters can be computed perturbatively 
because short-distance physics is perturbative. The above three-dimensional theory is super- 
renormalizable and has UV divergences associated with the parameters /is and /s. The 
coefficients of the other terms, however, have a simple finite relationship to the parameters of 
the original theory. The parameter /s represents the contributions of the non-zero Matsubara 



frequency modes to the free energy density (along with any associated UV counterterms of 
the three-dimensional theory). 

It is often conventional to rescale the field of the effective three-dimensional theory ( |1.12| ) 

1/2 



(TflT 
^1 (</>!+ ^02) (1.13) 



and write 



& = / (fx 



l|V0r + ^ 0^ + ^(0^)^ + ^3' 



:i.i4) 



2 I ^1 2 ^ 4! 

where = (0i, ^2) is a real 2-vector, 0^ = 0^ + 02, and 

'^bare = ^ , U= —^, J^^ = Pf^. (1.15) 

This is 0(2) scalar field theory in three dimensions. This form makes it easy to under- 
stand the scale at which physics becomes non-perturbative. The only parameters of the 
0-dependent part of ^3 are rbare and u. For fixed u, imagine finding the rbare(w) that cor- 
responds to the phase transition. Then all correlations at the phase transition can be con- 
sidered as determined by u. By dimensional analysis, the distance scale of non-perturbative 
physics is therefore 1/u. Z^ and Z^ turn out to be perturbatively close to 1, and so this 
scale is 1/m ~ A^/a by (|1.15| ). 

Perturbation theory in the three-dimensional theory ( p..l4[ ) is an expansion in u. Consider 
the dimensionless cost of each order of perturbation theory. By dimensional analysis, the 
contribution to that cost by physics at a momentum scale of order p must be order w/p. This 
means that perturbation theory breaks down for distance scales p~^ ^ u~^ . It also means 
that, thanks to ( |1 . 7] ) , perturbation theory works fine for distance scales p^^ ^ A at the 
transition, which are the distance scales of the non-zero Matsubara frequency modes. This 
is the reason that those modes can be treated perturbatively and a perturbative matching 
calculation is possible. 

To compute the number density at a given temperature and chemical potential, it is 
convenient to rewrite n = {ip*ip) as 

n = iPV)-'^ \nZ = -ipV)-' /^\ . (1.16) 



In the equivalent three-dimensional description ( |1.12|) , this becomes 



n ^ -i(3V)-" ^^' 




7^ (|V^o| ) + -TT- mw - -T^ ((^0^0) ) - ^^ , (1-17) 



9/i 2m djj, dfi m dfi 



^ The rescaling used here differs from that of Ref. |g] by a factor of Z,^ . This difference of convention is 
actuaUy moot since Z^ = 1 at second order. 



where all the expectations are taken in the purely three-dimensional theory (and UV- 
regularized, as necessary). The only approximation made here is ignoring corrections that 
would appear as higher- dimensional interactions in ( |1.12| ) which, as we've already said we 



will review later, do not affect the calculation of T^ at second order. Also, we shall see that 
the derivative (|V0p) term and the quartic ((V^oV^o)^) term are not relevant at second order, 
so that one may simply take 

A very similarly structured matching calculation was undertaken in Ref. 0] for matching 
the continuum 3-dimensional theory to a lattice 3-dimensional theory, where the purpose 
was to match the two theories order by order in the lattice spacing, in order to improve 
approach to the continuum limit in numerical simulations. The structure of that calculation, 
and the topology of the required perturbative diagrams, is identical to what we will need 
for the present task. The only substantial difference is that we will be evaluating those 
diagrams in the 3-1-1 dimensional theory, rather than 3-dimensional lattice theory, and that 
the perturbative expansion will correspond to an expansion in the scattering length (or an^^^) 
rather than the lattice spacing. There is additionally a trivial difference in presentation: Ref. 
Q did not keep track of 0-independent terms analogous to /s, which we use in ( |1.18|) , but 
instead discussed the matching of ip*4' directly. 

The generic technology of matching calculations was reviewed in the context of our 
current problem in section IV. A of Ref. 0. The idea is to perturbatively calculate an 
identical finite set of physical infrared quantities in the 3-1-1 dimensional and 3 dimensional 
theories, and then equate the answers to determine the parameters of the 3 dimensional 
effective theory. Each of the calculations must be IR-regulated, but the dependence on the 
choice of IR regulator will disappear in the final result of the matching (provided the same 
IR regulator is used for both theories). In our case, this is because matching is accounting 
for the differences of the two theories at the short distance scales (^ A) associated with the 
non-zero Matsubara frequency modes. For the specific purpose of the matching calculation, 
/x can be formally treated as a perturbation, along with the quartic interaction proportional 
to a. That is because the short distance scales ^ A are associated with energies per particle 
^ T, which is large compared to the chemical potential /i at (or very near) the transition. 
With /i treated perturbatively, the imaginary time Feynman rules are given in Table |. We 
will use the notation ko, Iq, po, ... to designate the Matsubara (imaginary time) frequencies 
associated with propagators with momenta k, 1, p, ..., and have introduced the short-hand 
notation 

^k^—. (1.19) 

In order to streamline calculations later on, it is useful to review the fact that the critical 
value fic{T) of chemical potential has a somewhat different dependence on non-perturbative 
physics than the critical value nc{T) of the density. Whereas nc{T) becomes non-perturbative 
at first order in interactions, non-perturbative effects do not enter the calculation of fic{T) 
until second order. (See Ref. for a discussion.) To determine f^ciT) to first order, it is 
adequate to do a purely perturbative calculation directly in the original 3-1-1 dimensional 
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3+1 dim. theory of ip 


3 dim. theory of V'o 


N. 


1 


2*' 

Svra 
m 


> » > 

X 


+)" 

87ra 
m 



TABLE I. Feynman rule s, appropriate for a matching calculation in a uni form B ose gas, for the original 
3+1 dimensional theory ( |l.9| ) of ^ and the effective three-dimensional theory ( 1.12[ ) of tpQ. The variable ko 
represents the Matsubara frequency of the field, while ojk = k'^ /2m. At finite temperature, loop frequencies 
Iq are summed over the discrete values Iq = 2iTnT with n any integer. In dimensional rcgularization with the 

MS renormalization scheme, a factor of M' — {e^^^'^M /^/A^Y should also be associated with each 4-point 
vertex but has not been explicitly shown above. 



= x' = 



-• — 



+ 




+ Oia^) 



FIG. 1. The vanishing of the inverse susceptibility at the phase transition, expressed in terms of 
diagrams at first order in a. 

theory. In particular, one can calculate the inverse susceptibility x~^ of 4' ^ind set it to 
zero to determine the transition point, as in Fig. |l|. Such a purely perturbative calculation 
would be inadequate at second order, for which one can instead marry perturbative matching 
calculations with non-perturbative results from the 3 dimensional theory p. 



II. UV REGULARIZATION 



Before starting a detailed calculation, we have to choose our convention for regulating 
ultraviolet divergences in our 3+1 and 3 dimensional effective theories. We will use di- 
mensional rcgularization, replacing the 3 spatial dimensions hj d = 3 — e dimensions. One 
convenience of this choice is that, in the 3+1 dimensional theory, loop corrections to the 
zero-energy 2—^2 scattering amplitude vanish at zero temperature and zero chemical poten- 
tial, so that there are no corrections to the identification of the a in ( p..9|) with the scattering 



length 15 



We will define UV-renormalized parameters using the modified minimum subtraction 
(MS) scheme, and we shall call the associated renormalization momentum scale M. At 
second order, this will not require any UV subtractions for the 3+1 dimensional theory, 
which can then be taken to be 



5t 



P 



dr I d? ^x 



r{dr 



2m J m 



where 



p7E/2 



A-K 



(2.1) 



(2.2) 



[The factor of e'^l"^ l\f^ in ( |2.2| ) is what distinguishes modified minimal subtraction (MS) 
from unmodified minimal subtraction (MS); the difference between the two schemes amounts 
to nothing more than a multiplicative redefinition of the renormalization scale.] The three- 



dimensional theory ( 1.12 ) is 



S3 = P d 



or equivalently 



S^ 



j3— e 



X 



^° ^ 2m 



/i3 ) ^0 + M^Z,^ (^o>o)' + M-73 
' m 



d^-'x 



l|V0r + ^0^ + M^^ 



l2\2 



+ M-^ J-g 



(2.3) 



(2.4) 



This theory is super-renormalizable and requires only a finite number of UV counter-terms. 
In particular, in the MS renormalization scheme, only the coefficient of 0^ is explicitly 
renormalized, with the exact relation 



_ 1 fu 



(2.5) 



between the bare coupling rbare and the renormalized coupling rj^^^M). 

The reader may wonder why we bother with the continuum 3-dimensional theory, 
since one will use a lattice-regulated 3-dimensional theory for actual computations of non- 
perturbative results. One could instead skip the continuum 3-dimensional theory and di- 
rectly match the 3+1 dimensional theory to the particular lattice theory used for a particular 
simulation. However, it is more convenient to split this matching into two steps: (1) 3+1 
dimensions to continuum 3 dimensions, and (2) continuum 3 dimensions to lattice 3 dimen- 
sions. The first step has the virtue of not depending on the details of how the theory is put 
on the lattice. 

III. THE MATCHING CALCULATION 

A. What we need 



There are two lattice simulation results of the three-dimensional 0(2) theory ( |1.14D that 
will turn out to be relevant to our evaluation of the number density n via (|1.18|) . Quoting 
values from Ref. 0], they are 



K 



A(0^)c 



u 



-0.001198(17), 



(3.1a) 



n = 



rf^M = u/3) 



u^ 



0.001920(2), 



(3.1b) 



where 



A(0^)c= (0^) 



\ 



(3.2) 



is the difference between the effective theory value of (0^), at the critical point, for the cases 
of (i) u small and (ii) the ideal gas m=O.0 Unlike (0^)c, the difference A(0^)c is an infrared 
quantity, independent of how the effective theory ( |1.14| ) is regularized in the ultraviolet. 
(Ref. 0] gives an independent and statistically compatible value of n but did not analyze 
Tc.) The reason that k and 7?. are pure numbers is dimensional analysis. If one picks the 
renormalization scale to be of order u, then, at the transition, the only parameter of the 
0(2) theory is the dimensionful parameter u. The dependence of A(0^)c and r^ on u is then 
determined by their dimensions. 

In dimensional regularization in the 3 dimensional theory, A(0^)c is the same as (0^)c- 
This is because, for the case u=0, the transition takes place at /is = 0, and then 



'h 



2 / -^ = Q 



(3.3) 



in dimensional regularization. (The last integral vanishes by dimensional analysis, since in 
dimensional regularization there is no dimensionful parameter to make up the dimensions of 
the integral.) Therefore, in the formula ( |1.18|) for n{T), we can replace the three-dimensional 

{"ipotpo) = ^^^"^"^{(f^) at the phase transition by Z^^mTnu, to obtain 



nc{T) 






(3.4) 



-Im=/^c(t) 



where we have used ( |1.15| ) for u. To evaluate nc{T) to second order in a oc u, we therefore 
need 



to first order in a and 



Za, Z^, and 









Mc 



Mc 



to second order. In evaluating the last quantity, we will find that we will also want fi^ to 
second order in a, which was computed in Ref. H. 



^ Technically, for the non-interacting case we must take the limit as the critical point is approached from 
negative /z. 




FIG. 2. The self energy at first order. 



B. Matching Z^ 

Z^ can be matched at first order by matching the infrared momentum dependence of the 
inverse Green function for ^o between the 3+1 and 3 dimensional theories. The one-loop 
contribution to the inverse Green function, shown in Fig. ^ is momentum independent. So 



0(a^ 



(3.5) 



where 0(a^) indicates corrections that are formally second order in perturbation theory. A 
slightly more detailed discussion is given in Ref. 0. 

In this paper, we will write 0(- ■ ■) when displaying the full parameter dependence of 
a correction (except possibly for logarithmic factors) and write 0(- • ■) when just showing 
the dependence on a particular parameter. So 32a^/A^ = 0(a^/A^) = 0[a^\ In matching 
calculations, where we are formally doing perturbation theory with IR regularization, (9(a") 
will just mean n-th order in perturbation theory. 

We now return to why we could drop the gradient term 



dZ^ 1 
(9/i 2m 



iV^oH 



dZ^ T 



dfi 2Z, 



i> 



(iv0r) 



(3.6) 



from the original matching formula ( |1.17| ) for n. At the transition, the dimensionally regu- 
lated (|V0p) must be 0{u'^) by dimensional analysis. This means that the expression ( p.6|) 
is at least third order in the interaction strength a and so irrelevant to our second-order 
calculation of ric- In fact, it is fifth order, since the dominant contribution to dZ^/dfi is 



C. Matching Za 

We will determine Za by perturbatively matching the four-point Green function of ipo at 
zero momentum. Specifically, to match Za at first order, we will compute the (amputated) 
diagrams of Fig. ^ in both the 3+1 and 3 dimensional theories. We need to choose an IR 
regulator for these two computations. For this case, we find the most convenient choice to 
be dimensional regularization, which will now be used to regulate both IR and UV infinities. 

In the 3+1 dimensional theory, the diagrams of Fig. |^ with zero external momenta give 
-MT(^) where 



r(4) 



ijca 



m 



Svra 



m 



2 Tp {ipo + ujp){-ipo + LUp) Tp {ipo + ujp 



C(a3). (3.7) 
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FIG. 3. The effective four-point interaction up to one loop. 



We introduce the short-hand notations 



Ti: 



po 



M' 



(27r)« 



(3.8) 
(3.9) 



where d=3—e is the number of spatial dimensions. The 3-dimensional theory result is the 
same as (p.7|) but with only zero-mode contributions included and factors of Za and Z^ 
inserted, so that 



r(4) 



-Za +i-Za ^ 



m 



m 



) 



T 



2Jp ui 



2T 



p UJ. 



p J 



+ 0{a^). 



(3.10) 



Equating the 3-1-1 dimensional result ( p.7| ) and 3 dimensional result ( p.lO| ), and keeping in 
mind that Za = I + 0{a) and Z^ = 1 + O(a^), we can now solve for Za through second 
order: 



1- 



Sna 



m 



1 — 5pQ,0 



+ 2 



2 Tp {ipo + ujp){-ipo + ujp) Tp {ipo + Up) 



1-5, 



Po,0 



+ 0{a'), 



(3.11) 



where 6ij is the Kronecker delta function. This expression is IR convergent, as it should be. 
In Appendix 0, we derive the dimensionally regulated results 



1 — Spofi 



p {ipo + ujp){-ipo + ujp) 



M'(3X 



-d 



c(i-i) 



(3.12) 



1-5, 



Po,0 



M^/?A-'^C(f-l) 



fp {ipo + UpY 
in d = 3 — e spatial dimensions. The result, after taking d ^ 3, is 

^.-l-l2C(l)^o($). 



(3.13) 



(3.14) 
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FIG. 4. Example of a ^-dependent contribution to the inverse susceptibility. 

For future reference, it's worth briefly stepping through the same calculation if we had 
separately evaluated the 3+1 dimensional and 3 dimensional contributions to the matching 
(the 1 and ^p^^o pieces of 1 — (5po,o), which are individually ill-defined without specifying a 
consistent IR regulator. We will often find it convenient to use dimensional regularization 
to regulate the IR (as well as the UV). The three-dimensional integrals are very simple in 
dimensional regularization, 

1 r d'^p 

-,ocy-f = 0, (3.15) 

p Up J p 

which follows by dimensional analysis. The full 3+1 dimensional piece would then be the 
same as the 1 — Sp^^Q results above. For example. 



i 



M^/5A-'^C(f-l)- (3.16) 



fp {ipQ + ujpY 

We are now in a position to explain why we could drop the quartic term 

dZa 27ia 



dfi m 



{{ro^oY) (3.17) 



from the original matching formula (|1.17 ) for n. From ( |3.14 ), or the diagrams of Fig. ^, we 



see that Za is /i-independent at first order. The leading /i-dependent contribution comes 
from graphs such as Fig. ^, which will produce an 0{fia) contribution to Za and so an 
0{a) contribution to dZa/dfi. There is an explicit a in (|3.17|) , which brings us up to 0{a^). 
Finally, there is the factor of 

((^o^o)^) = Z;^'m'T'{<p'). (3.18) 

At the transition, the dimensionally-regulated result for (0^) must be order -u^ by dimensional 



analysis. This means that the contribution (|3.17|) to n is fourth-order in the interaction 



strength and so irrelevant to a second-order calculation of nc(T). This argument is almost 
identical to a similar argument given in Ref. |^ for the matching of the 3 dimensional 
continuum theory to a 3 dimensional lattice theory. 

D. Matching the fi dependence of /13 

The matching of fi^ can be accomplished by computing the inverse susceptibility x~^ = 
r(2) [t^ both theories. At second order in a, this corresponds to the diagrams of Fig. |^. 
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FIG. 5. The inverse susceptibility through two-loop order. 

For the purpose of counting orders of a, we treat the chemical potential /i as 0{a). That's 
because we are ultimately interested in using the 3 dimensional effective theory at the phase 
transition - that is, for fi = fi^- The ideal gas result is /ic = 0, and the effect of interactions 
is that fie = 0{a). 

We discussed earlier that we need dfi^/dfi to first order. For this computation, we will 
find we can ignore the //-independent diagrams of Fig. ||. Specifically, Fig. |] gives 

- r(2) = 1^- -^IJ-i 7 7TT + (/i-independent) + 0{a^) (3.19) 

m Tp [ipQ + ujpY 

for the 3+1 dimensional theory and 

- r(2) = ;,3 _ Za^^ ^^T f ^ + (/i3-independent) + 0{a') (3.20) 

m Jp uo^ 

for the 3 dimensional theory. Equating the two results, order by order in a, yields 

yUs = Z^n + (yU-independent) + O{o?), (3-21) 



with 



Z, = l 1 ,. '"' 3.22 

m Tp [tpQ + Up)'' 



We then have 



^ = Z, + 0{a'). (3.23) 

The integral in p.22|) is the same as one of those encountered matching Za, and the result 



IS 



dn 



l-4C(i)^ + o(^). (3.24) 



A full result for ij,^ at exactly fi = fic was derived in Ref. 0. We will rederive the crux 
of that result in section |V| , since we will need it in a different analytic form than presented 
in the original derivation. 
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^o 




FIG. 6. Examples of diagrams which contribute to the free energy at up to third order. 







FIG. 7. Diagrams contributing to the /i-derivative of the free energy, through second order at the 
transition. Not shown is the tree-level contribution, dfs/dfi for the 3-dimensional theory. 

E. Matching the /x dependence of /a 

We now come to the more intricate part of our calculation, which is the calculation 
of dfs/dfi at the transition. We can match fs by computing the free energy density / in 
both theories. Some examples of diagrams which contribute to the /i-dependence of the free 
energy are shown in Fig. |^. We will find it convenient to instead compute the derivative 
directly, a diagrammatic representation of which is given by Figs. 0-§, which contain all 
diagrams contributing through O(a^) in perturbation theory. The filled circles are still 
associated with a factor of /z, but the crosses, which represent factors of /i that have been 
hit by a derivative, are associated with a factor of 1. One can also think of these diagrams as 
representing mixing of the operator 0^ (represented by the crosses) with the unit operator, 
which were the words used to describe analogous calculations in Ref. [ffl. 





FIG. 8. Contributions to the /^-derivative of the free energy which cancel at second order at the 
transition. 
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FIG. 9. The remaining second-order contributions to the /z-derivative of the free energy, which also 
cancel at second order at the transition. 

At the transition, the diagrams of Fig. ^ cancel at second order in a, as do the diagrams 
of Fig. 1^, because of the first-order relation of Fig. |I|. So we may simplify our task by focusing 
on just the diagrams of Fig. 0. 

At this many loops, organizing all the terms of the matching calculation becomes tedious 
unless one from the start uses dimensional regularization for the IR (as well as the UV). In 
the 3-dimensional theory, the loop integrals in Fig. |^ then all vanish by dimensional analysis, 
leaving us with the formal result 



dfjL 



dh 
djj, 



(3.25) 



for dimensionally regulated perturbation theory. In the 3+1 dimensional theory, we have 

1 . MK^ 1 






Mc 



p {ipo + ujp) 



at the transition, where 



/^c 



(/^c 



Sna 



/^i^': 



p {ipo + ujpY 



+ (^ribbaii + C(a^ 



m 



Tp {ipo + 



Ur, 



is the first-order result for /ic, determined by Fig. |I], and 

2 . 1 



Sn^h: 



■all 



5 Tray 



m 



J TpQK {ipo + ujpy{iqo + ujg){iko + uJk){i{po + qo + ko) + a;. 



p+q+kj 



(3.26) 



(3.27) 



(3.28) 



is the contribution of the last diagram of Fig. ^, which we refer to as the basketball diagram.^] 
In section |V], we will show how to evaluate this integral in dimensional regularization, with 
the result that 



6n 



bball 



327ra2 



2 -^2 + (2i^i + ln2)C(|)-i^ 



c(i; 



+ 0{e) 



(3.29) 



where Ki and K2 are numerical constants defined byQ 



^ Historical note: The origin of this terminology in the literature is associated with the physical appearance 
of old American Basketball Association basketballs, not current National Basketball Association ones. 



'^ For efficient numerical evaluation of Ki, it's useful to make a change of integration variables in ( |3.30 ) 
such as u = \/l ^ tj which makes the integrand analytic at both endpoints. 
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At: Jo t 



Lii/2(t) 



-Kt 



_ 1 /■! ds dt 



l-i 
iLii/2(s)Lii/2(t)LLi/2(st) 



0.630 568 207 496 069, 



(3.30) 



sty/n 



2(1 - st)3/2 

-0.135 083 353 73. 



vr 



1-s 



c(i; 



vr 



1-t 



+ C(|) 



Li,y(2;) is the polylogarithm function defined by 



n=l 



n" 



(3.31) 



(3.32) 



[The polylogarithm function is often called gy{z) in statistical mechanics.] 
The other integrals in the matching ( p.26| ) are given by (|3.16| ) and 



M'\-'^C{i) 



(3.33) 



IP Wo + ^p 
which is discussed in Appendix 0. Putting it all together and taking d ^ 3, 

- (f^) = C(i) A-' + C(i) A-^ (/3/ic - /3/iW) + <5nbbaii + 0{a') . (3.34) 

In Ref. 0, a second-order matching calculation was carried out for ft^ that determined /ic 
in terms of the critical value of r in the 3 dimensional theory. The result was 



/3/ic = /?/iW + 



32™^ 



In(MA) + Ci- 727r' 



rc(M) 



u^ 



(3.35) 



with Ci a numerical constant that was given in that reference in terms of a somewhat 



inelegant double integral. In section 0, we show that Ci can also be expressed as 



Ci = i - i ln(327r) - 2Ki. (3.36) 

Combining the last several formulas, and choosing M = u/3 to make contact with the quoted 
lattice measurement ( p.lb|) of Vc , we have 






^1 =C(|)A-^ + 



Mc 



32na^ 



X 



In - ) + i ln(1287r3) + i - 72n^n C(| 



vr 



Ko 



In 2 



ai: 



-hy 



(3.37) 
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IV. FINAL RESULTS 



We now have all the elements we need for nc(T) as determined by 



n, 



.(T) = X-^ [bo + biaX'^ + [b'^ ln{aX-^) + b'^]a^X~^ + 0{a^X-^)} 



The result is 



(4.1) 



with 



60 = C(|), 

61 = 1927r^K, 

6'2 = 327rC(i), 

62' = 327r| [i ln(1287r3) + i - 727^27^ - gevr^/tl C(| 



K. 



In 2 



2 "' 2^/? 
Inverting this formula, 

T,{n) = To{n) 1 1 + Cian^^^ + [4 \n{an^/^) + 4'] a^n^^^ + 0{a^n) 

with 

„ _ 2L-4/3, 
Ci — —3 Oq Oi , 

'-2 ~ 3 "0 "2 ' 

// _ 2 /,-5/3,// 7 L-8/3, 2 1 ^' It, /, 
^2 — ~ 3 Oq O2 + g Oq f ;^ — 3 C2 m Oo , 



C(^ 



(4.2a) 
(4.2b) 
(4.2c) 

(4.2d) 



(4.3) 



(4.4) 
(4.5) 
(4.6) 



where the ln(6o) term arises because we have changed the argument of the log from ln(a/A) 
to ln(an^/^). Putting in the lattice results of ( |3.1|) , we get the numerical values 



61 = -7.1(1), 62 = -146.8108- 



b'^ = -587(2), 



(4.7) 



ci = 1.32(2), 



4 



19.7518- 



4' 



75.7(4) 



(4.8) 



These are our final results. 

In our discussion of effective theories, we left out a variety of corrections, such as 3-body 
interactions or cross-section energy dependence in the original 3+1 dimensional theory, or 
(j)^ and higher-dimensional operators in the 3-dimensional effective theory. Ref. ^ contains 
a detailed discussion of the parametric size of the resulting corrections to Tc{Np) for an 
arbitrarily wide harmonic trap. The same analysis holds for Tc{n), with the result that 
there are no corrections at second order. However, a third-order result would depend not 
only on the scattering length a but also on the effective range of the two-body scattering 
potential. 

The relative size of the second-order result obviously depends on the diluteness of the 
gas and the value of the scattering length, which will vary from experiment to experiment. 
However, just for fun, let us produce numbers for one particular case. In 1996, Ensher et 
al. Il^] studied the BEC transition for for dilute gases of ®^Rb atoms in the F=2 hyperfine 
state, trapped in a harmonic trap. The relevant scattering length is a = (103 ± 5) a^ ||17|| , 



17 



where oq = 0.0529177 nm is the Bohr radius. Their transition was at T ~ 280 nK. These 
parameters correspond to a/Ao — 0.015. Let's now consider a theoretical prediction of the 
central density n (which, sadly, is not directly accessible experimentally). Eq. ( |4.1| ) gives 
a first-order correction to the ideal gas result no{T) of roughly —4.2% and a second-order 
correction of roughly +0.2%. As discussed in Ref. [0, where a similar analysis is made of 
Tc{Np), this particular trap may not actually be wide enough for our second-order result to 
be valid. 



V. EVALUATING ^ribbaii WITH DIMENSIONAL REGULATION 



We will now explain how to obtain the result ( |3.29| ) for the basketball diagram in dimen- 



sional regularization. We do not know how to carry out the integrations for the basketball 
( p.28| ) in arbitrary dimension d. Instead, we will find it convenient to manipulate the inte- 
grals into a form where we can dispense with dimensional regularization and do integrations 
in three dimensions. We are currently relying on dimensional regularization to regulate the 
IR and UV divergences of our calculation. Our first step on the road to dispensing with 
dimensional regularization of the basketball diagram will be to convert the IR role of di- 
mensional regularization to mass regularization, which will turn out to be more convenient 
for this particular calculation. To make this conversion consistently, we will first relate the 
basketball diagram to a combination of diagrams that is IR convergent. In particular, con- 
sider the combination of diagrams Sn'y^^^Yi depicted in Figs, [l^ and ^ (ignoring all Z^ and 
Za factors). In equations, this is 

^ / 1 / Svra 

"^bball 



2 V m 



where 



PQK {ipo + UJpYiiqo + UJq){ikQ + t^fc)(«(PO + go + ^o) + tUp+q+k) 

-T=^ / —. ^ 1 - Ansu„(0) T [ \, (5.1) 



Ans,n(o) = - ' ^ 



2 V m 



QK {iqo + ujq)iiko + uJk){,i{,qQ + h) + ^q+k) 
1 



rpZ 



qk LdqiOk^d^]:^. 



(5.2) 



is the difference between the "sunset" diagram contributions to the self-energy at zero mo- 
mentum in the 3+1 dimensional and 3 dimensional theories. The above can be combined 
into 

r / _^ ( 87^a^^ ^ 1 — 5po,o^<jo,o^fco,o 



'bball 



2 V m J TpQK {ipo + ujpy{iqQ + ujg){iko + Uk) 

(5.3) 



1 ^Po,0 



.i{po + qo + ko) + Wp+q+k iiqo + ^o) + '^q+k- 

This expression is IR convergent. (Actually, it is absolutely convergent in the IR only if one 
first averages the integrand over p -^ —p. This is a technical point that won't have any 
impact, and we shall implicitly assume such averaging wherever required.) 



3+1 dim. 



3 dim. 



3 dim. 



5n^^>= 





+ 




FIG. 10. IR convergent combination of diagrams corresponding to (Jribb, 



all- 



3+1 dim. 



3 dim. 




p=0 + 




P=o = -An,„„(0) 



FIG. 11. A piece of the second-order correction to the relationship between /i and r, which is to be 
understood to be evaluated at zero external momentum. 



In fact, in dimensional regularization, the original basketball ( 3.2S ) is equal to the IR- 
convergent version ( p.l| ). That's because the explicit three-dimensional integrals of the 
terms that have been added to ^ntbaii in (|5.1| ) all vanish in dimensional regularization by 
dimensional analysis, as in ( |3.15| ). 

So we can now focus on ( pl3|) . Because it is IR convergent, it will not be changed if we 
add an infinitesimal negative chemical potential 5/i to the propagators. We then re-expand 
the terms into the form of ( p.l| ), where the IR divergences of the individual terms are now 
regulated by 5ji: 



(^^bbaii = lim 



Eq. (5.1) with uj — > uj^ and 11 -^ 11 



(5.4) 



An+^(0) = [Eq. (5.2) with u ^ u 
where 

U^ = OJ — 5jJL. 

Henceforth, the limit (5/x -^ 0^ will be implicit when we discuss (5nbbaii- 



(5.5) 



(5.6) 



A. The mass-regulated basketball 



We will now focus on the first term of (5.41 



5ni = 



ina 



m 



+ \' 



PQK {ipo + uj+y{iqo + uj+){iko + cu^)(i(po + Qo + h) + ^p+q+k 



(5.7) 
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which is just our original basketball diagram with an infinitesimal chemical potential 6^, 
and with the UV still regulated with dimensional regularization. This can be expressed as 
the derivative with respect to 5/i of the corresponding contribution 



Pi = 7; 



BvraN 
m ) TpQK {ipo + uj+){iqQ + uJ+){ikQ + uj^){i{pQ + go + ^^o) + ^p 



(5.8) 



+q+kJ 



to the pressure, corresponding to the last diagram of Fig. |^. For arbitrary chemical potential 
(not just the infinitesimal case), this contribution was derived in Ref. ^ for UV dimensional 
regularization, and simply reproduces a corresponding portion of a 1957 calculation by 
Huang, Yang, and Luttinger |jT8|, which used a different UV regulator. The result (after 
taking d — > 3) is 



z 



i+j+k 



Qrp 2 00 00 00 

^' ^ ^^§§£F+ A:)(J + A:)(uA:)V2' 
where z is the corresponding fugacity. In our case. 



(5.9) 



(5.10) 



which is infinitesimally less than one and is serving the role of an IR regulator. Differentiation 

gives 



with 



Sni = -^S{z) 






(5.11) 



(5.12) 



Our task is to find the 5/i — i> (^ ^ 1 ) behavior of this sum, extracting any divergences 
and the finite remainder. It will be convenient to define 



a 



-pSfx, 



so that the limit of interest is a ^ 
Now rewrite the sum as 



ijk 



+ j + k)z'+^+'' /-i ds dt ■ 



{ijkyn 



+k.j+k 



S t 



s'+H^ 



fdsdt^ ^i+k^j+k^i+j+k 



1/2 



1/2- 



J^. 



+ 



.^3 . 



If one does the sums before the integrals, the sums can then be factorized, giving 

"i ds dt r 



S{z) 



s t 



(5.13) 



(5.14) 



2 Li_i/2{sz) Lii/2{tz) Lii/2{stz) + Lii/2(sz) Lii/2(tz) Li_i/2(stz) . (5.15) 
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If we removed the mass regularization by simply setting z = 1, the above integral would 
have singularities associated both with (i) s — *> 1 with t fixed, and (ii) s and t approaching 1 
simultaneously. The first type of divergence arises only from the first term in the regulated 
integral (|5.15|) and can be eliminated by integrating this term by parts. Using 



— U^{sz) = -U^_i{sz), 
as s 



integration by parts gives 



"^ds 



/■i (ig fi dg 

/ — Li_i/2(sz) Ui/2{stz) = Lii/2(2:) Lii/2(te) - / — Lii/2(s2;) U_i/2{stz). 

Application to ( |5.15| ) then yields 

rl dg dt 
S{z) = Si{z) - — — Ui/2{sz) Lii/2(t2;) U_i/2{stz), 
Jo s t 



where 



r^ dt [ 1 2 

Si{z) = 2Uii2{z) — Lii/2(t2;' 
Jo t ^ 



(5.16) 



(5.17) 



(5.18) 



(5.19) 



To analyze the behavior of this expression as 2; ^ 1, it will be useful to have the following 

(5.20) 



n=0 



series expansion of polylogarithms |19], 

00 
Li,(e-") = r(l-z/)a''-i + 5] 

The relevant special cases are 

Li+i/2(z) = 



n\ 



■C(^-n)«^ 



TT 



Li-l/2(2) 



I- Z 



+ a\)+o{VT 



+ 0{{l-z)-^") 



- + C(|) + 0(v^), 



+ 0(a-i/2). 



TT 



2(1-^)3/2 ' vv J ' 2a3/2 

Let's separate out the regularization dependence of (|5.19|) for Si by writing 



(5.21) 
(5.22) 



^dt T. , , n 2 

— Lii/2(te 

t L 



1 -adt 



+ 



^dt 



Lii/2(tz) 



vrt 



1 — -Zt JO t 

ln(l-2;)+ / - Lii/2(tz) 

z Jo t ' 



1- zt 
2 nt 



1 - zt 



(5.23) 



If we just wanted an expression through 0((1 — 2;)°) = 0{q^), we could now set z = 1 in the 
very last integral. However, the Lii/2(2;) factor in (|5.19|) has an 0(a~^/^) singularity, which 
means we need the expansion of ( p.23[ ) through 0(a^/^). This term is easily obtained by 
differentiating the last integral in ( ^.23| ) with respect to z and then analyzing the dominant 
piece of the result, which is a singularity (cut off by 2) as t ^ 1. The result is 



^ dt V. , , n 2 

— Lli/2(tZ 
t L 



-TT In a + 4777^1 - 4C(i) Vvra + 0(a), 
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(5.24) 



where Ki is as defined in ( |3.30|) . So 



Si 



- + C(i) (-27r Ina + SrcKi) - SnCH) + O(v^). 



(5.25) 



Now let's return to the expression (|5.18| ) for the sum S and work on isolating the re- 
maining (regulated) divergences associated with s and t simultaneously approaching 1. We 
isolate the singular pieces of the integrand in ( |5.18|) by writing 



S{z) = Si{z) + S2{z) - 4itK2 + O(v^), 



(5.26) 



Soiz) = 



1 

(Xo Lilt' „ , > ,-, /", 

2(1 - stz)3/2 



n 



1 — sz 



+ C(i) 



TT 



1-tz 



+ak: 



(5.27) 



with K2 defined as in ( p.31|) . [The C(2)^ term in 5*2 doesn't actually give a singular piece, but 
including it makes the remainder (|3.31| ) a little more compact to write.] Explicit integration 
gives 



So 



7]-3/2 ]y. A r n 2 

- + 271 Ca) ln(4a) - V^ C(i) ln4 + 271^/2 + O(v^), 



(5.28) 



and so 



Si + S2 = 47r ,/ i ln(2a) + 2Ki 



+ 



+ {2K,~2 + ln2)Ca: 



In 2 



ai: 



O(v^). (5.29) 



The combination of ( [5. Ill ), ( p.2(j| ) and ( |5.29| ) then gives our result for the mass-regulated 
3+1 dimensional basketball. 



6ni 



32™^ 
~1^ 



l\n{2a) + 2K^ 



In 2 



K2 + i2K,-2 + \n2)CC2)-^[a-2)\ ^^(v^)' (5-30) 



where K2 is defined as in ( |3.31|) . 



B. The remaining pieces 



The second diagram of Fig. 10, when mass regulated, gives a contribution to (Jribbaii of 



5n2 



m J /--'- ('/. .+^2, ,+,,+, 



pqk [UJ^j-U^UJi, ^p+q+k 



-2^^TT^a^m^T^ 



pqk (p2 + A^2)2(g2 + M^){k^ + M^){\p + q + k|2 + M^) ' 



(5.31) 



where 
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M"^ = -2171 5 IX 



Aira 
1^ 



(5.32) 



is infinitesimal. It is not strictly necessary to calculate this term because, by dimension 
analysis, the last integral is proportional to 1/A^ (in three dimensions) and so will only con- 
tribute to the cancellation of IR divergences in ( p.l| ) and not the finite remainder. However, 
it's reassuring to check the cancellation. From Ref. ||TT[|,p| 



1 



pqk (p2 + M^){q^ + M^)ik^ + -M2)(|p + q + k|2 + M^) 



IQn^ 



In^^ +ln2-2 

2e 2 V M / 



+ 0(e). 



(5.33) 



Now differentiate with respect to A^^ to obtain minus 4 times the corresponding integral in 
( p3l| ). Then 



6n2 



327ra2 Rf 



A5 



_l 3^Jl6^^ 



, .... - , , In 2 - - 

2e 4 VM^AV 2. 



(5.34) 



Finally, we need the last diagram of Fig. [1^ with mass regularization, corresponding to 

1 



6n3 = - 
In dimensional regularization, 

1 _r(2-i) 



-An+„(0)T 



HY 



(5.35) 



M 



d-i 



giving 



p(p2 + A<2)2 (4^)d/2 

An+,(o) /Tf 



%TlM 



eln 



2M 



M\ 



M J 



OU' 



(5.36) 



(5n,, 



A3T 



e / 167ra 
2 '^VM2A2 



)+0(6^ 



(5.37) 



Because of the a ^^^ in this equation, we will need the result for the mass- regulated An^j^(O) 
through 0{y/a.). The 0(a°) piece, which is independent of the IR regulator, was calculated 



in Ref. ||§] as part of calculating /Xc, and we will rederive it in section 0. The result is 

An3un(0) = -^^^ |1 + In(MA) + Cij , (5.38) 

with Ci as in ( |3.3(j| ). To get the 0{^/a) piece in the mass- regulated version, we start with 



the integrals corresponding to evaluating Fig. 11 with an infinitesimal chemical potential 



An+„(o) 



1 



57ra 



1 ~ ^qofl^kofi 



m J Tqk {iqo + uj+){iko + cj^)(i(go + ^o) + t^q+k) 



(5.39) 



See also Ref. |20|, which has a useful collection of dimensionally regulated three dimensional integrals. 
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Now differentiate with respect to a = — /3 5/i, 



d_ 
da 



An+,(0) = -T 



STraV 
m J 



1 ~ ^qofl^kofl 



QK {iqo + uJ+Yiiko + a;^)(i(go + h) + t^, 



q+kJ 



+ 



1^ 



^go,0^fco,0 



2fQK {iqo + uj+){iko + uj^){i{qo + ko) + uj, 



q+kJ 



. (5.40) 



We want to find the divergent 0(a~^/^) pieces of this expression in order to obtain the 
0{^/a) piece of ( |5.39| ). The first integral in ( p.4UD has an IR divergence associated with 
go = and g — > 0. In this limit, it can be simplified to 



T 



^ko,0 



{uJ+yTk (iko + ujl 



+ 0(«°). 



(5.41) 



The infinitesimal chemical potential can be dropped in the k integral, since fco 7^ cuts 
off the infrared, and then the integral is given by ( p.l3|) . The q integral is proportional to 
(^). So 



■'fco.O 



m'^Ci 



2 /•/'l-i 
2> 



ci {uj+y Tk {iko + uj'^y AnX'^y/Tca 



0{V^)+0{e 



(5.42) 



The second integral in (|5.39| ) can be evaluated similarly by first making the change of 
variables Q -^ Q — K, to get 



1 — Snn.oS 



qofi'Jkofi 



QK {iqo + uJ+Y{iko + u;^)(i(go - ko) + uj, 



q-kJ 



T 



LiuJ+yfK 



i-s, 



kofl 






(uj+yfK {iko + ujt){-ikQ + ut) 
2^ +0(«°) + 0(e). 



+ 0(a°) 



(5.43) 



Putting it all together to get 9An+jj(0)/9a and then integrating gives 

647ra^T ,,i, [a 



An+,(0) = An,,,(0) - ^-^^ C{\)xl- + 0{a) + 0(e). 

A 



(5.44) 



Combining with ( |5.37 ) and ( 5.38 ) 
327ra2 



Sn-^ 



A5 



- + ln(MA) + C,--l„(^ 



+ 2c(i: 



Combining ( ^.30| ), ( ^.34|) , and ( ^.45| ) for the three diagrams of Fig. |T0 

'^'^bbaii = 5ni + 5n2 + (5n3, 



as 



(5.45) 



(5.46) 



together with ( p.36|) for Ci , then yields our final result ( p.29|) for the dimensionally regulated 
basketball (Jribbaii- AH the IR divergences cancel, as they should. 
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VI. REDERIVATION OF fi^ AND AHsunCO) 

Second-order matching results for /ic and AHsunlO) are derived in Ref. |^. However, for 
the purposes of this paper, we want them expressed in terms of the same sorts of polylog- 
arithm integrals that we used in our evaluation of the basketball diagram. Here, we shall 
show how to obtain that form. 



Start with Fig. 11 and Eq. (|5.2|) for AH(0). This expression is infrared convergent and 



so independent of the choice of IR regulator. [We will not concern ourselves here with 



vanishing corrections, such as the 0{^/a) piece that was important in section [VB| .] As we 
have done before, it will be convenient to nonetheless introduce an IR regulator and evaluate 
the 3+1 dimensional and 3 dimensional pieces separately. For these particular diagrams, the 
convenient choice of IR regulator will be to introduce an infinitesimal chemical potential on 
just two of the three internal propagators: 

AH,,,(0) = -Hi^+i)(0) + H(i(0) 

^1/ 87ra\2^ 1 - ^90,0^0,0 ,g ^^ 

2V mj Tqk {iqo + ujg){iko + uj^){i{qo + ko)+uj:^^^)' 

The frequency sums in the 3+1 dimensional piece can be evaluated using standard contour 
tricks as in Ref. [pi to yield 



1 /87ra\2 f ,,,,,, /3c^+ _ g/3K+c^+) 

f n{ug) n{ujl) n{uj,^) — ^ 

qkl Ui 



-^ltl{^) = \r^) / n{u,)n{^t)n{^t)-^^ ^ (2.)^5(^)(1 - q - k). 



1 /Svra^ 2 



2 \ m J Jqki 



^(2.)^^M(l-q-k)^ (6.2) 



Note that the infinitesimal chemical potential (5/i cancels out in the denominator uj'^—ujq—uj^. 
This fact will simplify our analysis later on, and it is the reason that we chose to put regulator 
masses on only two of the internal lines rather than all three. Now apply a redundant 
principal part (P.P.) prescription to this denominator so that we can separately evaluate 
the integrals of each term. The first term in ( |6.2|) vanishes on angular integration. The last 
term in ( |6.2[ ), involving just one factor of n, is 0(e) in dimensional regularization for the 
same reasons discussed in Ref. 0, which are that 

pp(2.)^.(-'(l-q-k) ^ ^^^j ^^ 3j 

qk UJi — UJq— UJk 

and that the remaining 1 integration is convergent and cannot generate a compensating 1/e. 
Exchanging integration variables in the second term, we are then left with 

2 \ m J Jqkl L '^ J (^^ _ ^^ _ (^^ 

(6.4) 
25 



This integral has no UV divergences, and so we can set c? = 3 in the integral. Now expand 
the Bose distribution functions as a series in the regulator fugacity 2;, 



a=l 



^a^-aPup 



(6.5) 



giving 

-ni^rno) 



/ -1. v/qkl 



2 \ m 



ab 



^l — ^q — ^k 



Rescaling integration variables to be dimensionless, 



-n(^rno) 






^(^«+^ + /)/„, + 0(e), 



ab 



where 



jqki |(r — q^ — k'') 



(6.6) 



(6.7) 



(6.8) 



Using the methods of Appendix A of Ref. 0, one may evaluate this integral, obtaining 

1 



(a + h)\fah 
We can now extract the IR divergences of the sum using the same method as in section 



V^ 



ab 



^dt r 
T 



Lii/2(2;t)] = 7rln« -47ri^i + 0(v/a). 



^ /J,fe = - /' y Lii/2(t) Lii/2(^t) = vr In (^) - 47riri + O(v^), 



So 



- nSr^(O) = ^^ {ln(f ) - AK,] + OiV^) + 0{e) 
The corresponding result in the 3-dimensional theory is 



n(i(o) 



A2 



-J(-M), 



(6.10) 
(6.11) 

(6.12) 
(6.13) 



where 21 



J{M) ^ I 



1 , M 1 
2^ + ^"2A?+2 



/qkg2(fc2 + _/t<2)(|q + k|2 + A<2) (47r)2[2e ''"2M ' 2j^^^^^' ^*''^^'' 

Putting everything together, one obtains the previously quoted result (|5.38|) for Ansun(O), 
with the constant Ci given by ( [3.361 ). This new form of Ci, which is equal to the value 
derived in Ref. 0, can then also be used in the result of Ref. for yUc, which we quoted in 
(PI). 
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VII. THE SECOND-ORDER LOGARITHM FOR ARBITRARY N 

As mentioned earlier, Holzmann, Baym, and Laloe argued for the existence of the loga- 
rithmic term at second order. In order to make their general argument more concrete, they 
also presented an approximate large A^ calculation of the coefficient c'2 of that logarithm. 
It's interesting to compare exact results for the coefficient Cg to their approximate large A^ 
calculation. For this reason, let us consider generalizing our 3+1 dimensional theory ( |1.9|) 
to a theory with A'^c complex fields with U(A''c) symmetry: 



S^= f dT fd^x 



E*-(9-iv=-'')* + ^fev«.. 



(7.1) 



The corresponding 3 dimensional theory can be considered a theory of N=2Nc real fields and 
has 0{N) symmetry. The action is again ( [1.12 ), with u as before, but now 0^ = </>i + - ■ ■ + 0Ar- 



For A^c ~ 1; we reviewed before that the dimensionless cost of each order of perturbation 
theory gets a contribution of order u/p from physics at momentum scale p. For large A'c, the 
contribution is order Ncu/p. (See Refs. PJ2^ for a discussion of large A^.) In both cases, the 
momentum scale of non-perturbative physics can therefore be characterized as order N^u, 
and the condition for the theory to be perturbative at the scale p ~ A~^ associated with 
non-zero Matsubara modes is N^uX <^ 1. This generalizes the previous condition ( [1.7|) for 
useful expansions of Tc{n) and the applicability of perturbative matching to 

N^a < A (7.2) 

at the transition. We shall assume this in what follows. 

The second-order logarithm in our calculation of ndT) arose from the second diagram 
of Fig. 1^, which is proportional to fic at the transition. More specifically, it comes from the 
sunset diagram Ansun(O) contribution to /ic (|3.35|). (The sunset diagram is depicted in Fig. 



[TT].) The value of the sunset diagram for general A"c is simply the value for A"c = 1 multiplied 

by 

N^ + 1 N + 2 , , 

' = -^—- (7.3) 



So, ignoring non-logarithmic second order corrections, Eq. (|3.35|) for fic is modified to 



/JA. = /?/.!" + ^ 2?^ ln(A«) + .... (7.4) 

Now recall that for the case A'c = 1 analyzed in the rest of this paper, we chose the renor- 
malization momentum scale M to be of order the momentum scale u for non-perturbative 
physics. The reason goes back to section |II1 7\| : such a choice makes the critical value of 
rc(M) proportional to u"^ by dimensional analysis. That meant that the the rc{M)/u'^ term 
in the formula ( p.35[ ) for /ic really gives a straight C(a^) correction without any logarithmic 
enhancements. So, if we want to remove the possibility of implicit logarithms in the r^ 
term hiding in the ■ ■ ■ in ( |7.4| ), we should again choose M of order the momentum scale of 
non-perturbative physics, which in the present context is 
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M ~ N^u. (7.5) 

So 

^ m (AT, + 1) 32™^ /A/'.aX ^ /A^ea^X ^ ^ 

/3/ic = /3/xW + ^ \ ^ ^^ ^""VJ^) + ^ ( 3^ ) ' ^^-^^ 

where, for this section only, the notation 0{Nca^ /\^) is meant to assert that there are no 
additional factors of logarithms. 

The second diagram of Fig. |^ gives a factor of Nc as well as the factor of /Xc- The ideal 
gas result no(T), depicted by the first diagram, also has a factor of N^. The generalization 
of (|1D for n^{T) is 



n,{T) = N,X-^ I bo + haX-' + li!i±ll byX~''[\n{N,aX-^) + b'^] + OiN.a^X-^) 



(7.7) 



where bo, bi, and 62 are as in ([4.2|) but with 



The solution for Tc{n) can be written in the form 



A((/)')c . ^ 

K= l^ '\ (7.8) 



ATc 1287r'^fi; a ^^ , /A^cOX „ / a 



3, , ■ ^TflM^ +0T2h (7-9) 



To C(2) -^0 Aq V Ao / \A 

where 

32vrC(|; 

'3^ 



^ = -<^= + ''Tcf '""' 



and Ao = A(To). In the large A^ limit, k -^ — l/967r^ |J2^ and the coefficient A of the 
logarithm becomes 



327rC(| ) 
3C(|) 



^ -^ -^c ^^;l^ ^ 18.7327 A/'c. (7.11) 



This is the exact large A^ result for this coefficient if large A^ is defined to mean taking 
a large number of fields in the original 3+1 dimensional theory, as we have above. For 
comparison, the calculation of Ref. [^ was an approximate large A^ calculation made in the 
O(A^) three-dimensional theory using a rough physically- motivated UV momentum cut-off 
of A «i (27r)^/^/A. Their approximate result for A, expressed in terms of A^c = 2A^, was 

2567r 2567r 

Ab.,. ,5,, ^ A'e 3^j3y^ « ^ 3C(j)(2.)V^ " '"■''''■ '"'' 



This differs by roughly a factor of two from ( [r.ll| ). Amusingly, the approximation ( [7.12| ) of 



Ref. 1^ does well if naively applied to A'c = 1, giving A(Ref. |§) ~ 40.9. In contrast, the exact 
answer (|7.1CI|) derived in this paper is A ^ 37.4565. 
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APPENDIX A: SOME DIMENSIONALLY REGULATED INTEGRALS 

In this appendix, we show how to do basic single-momentum integrals used in the main 
text, such as ( |3.12| ), (|3.13|) , and (|3.33| ). We will do the po frequency sums first, using standard 
contour tricks. One can just as easily get the same results by doing the p integrations first 
(though one must be careful in that case about cuts). 

Start with the basic integral 



t 



"n{uj,-f,) + l]. (Al) 



ptpo + ujp- jJ. 

The +^ comes from the contour integration at infinity but can be ignored, since its integral 
vanishes in dimensional regularization. The Bose distribution function can then be expanded 
as in (|6.5|). If each term is integrated in d=3—e dimensions, using the t = case of 



d^p ,^_v _ ^-^^^*^/^r (^) 



pt^-Ap ^ VZ^ ^ ^^2) 



(27r)^^ (47r)<^/2r (f ) ' 



one obtains 



1 



ptpo + Up- fi 



M'X-^U.iz), (A3) 



where z = exp(/3yu). If we set /x = in dimensions where the result will still converge (and 
then analytically continue to other dimensions), we get the integral ( p.33|) used in the main 



text. If we first differentiate with respect to /i and then set /i = 0, we get ( |3.16| ) [which, as 
discussed in the text, is equivalent to ( p.l3| ) in dimensional regularization]. 
We can do the integral ( |3.12| ) of the main text as 



1 



p {ipo + ujp){-ipo + ujp) Tp 2ujp 



1 1 

+ 



n{ujp) + 1 _ f n{ujp) 



tpo + ujp -tpQ + UJp\ Jp Up Jp Up 

(A4) 



The last integral can also be done by expanding n as above, using the t= — 2 case of ( |A2D , 
with the result quoted in the text. 
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